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Wc describe the change of the spatial distribution of the state of polarisation occurring 
during two-dimensional imaging through a multilayer and in particular through a layered 
metallic flat lens. Linear or circular polarisation of incident light is not preserved due to 
the difference in the amplitude transfer functions for the TM and TE polarisations. In 
effect, the transfer function and the point spread function that characterize 2D imaging 
through a multilayer both have a matrix form and cross-polarisation coupling is observed 
for spatially modulated beams with a linear or circular incident polarisation. The point 
spread function in a matrix form is used to characterise the resolution of the superlens 
for different polarisation states. We demonstrate how the 2D PSF may be used to design 
a simple diffractive nanoelement consisting of two radial slits. The structure assures the 
separation of non-diffracting radial beams originating from two slits in the mask and 
exhibits an interesting property of a backward power flow in between the two rings, (c) 
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1. INTRODUCTION 

In 2000, Pendry [T] showed that a simple silver slab is capable of imaging in the near-field with resolution 
beyond the diffraction limit for the near-UV wavelength range. The operation of a similar but asymmetric [2] 
flat lens was later experimentally verified [31I3]- In these works, the mechanism leading to super-resolution 
is related to the appearance of a surface plasmon polariton (SPP) mode, which enables transfer of the 
evanescent part of the spatial spectrum through the slab. Other ways of achieving superlensing include the 
use of photonic crystals showing negative refraction [5H5] . Both losses and the cut-off wavelength of the SPP 
mode limit the operation range to sub- wavelength distances only [TU1[TT| . This obstacle has been tackled to 
some degree in several ways. One of them is by breaking the silver slab into a metallic multilayer [2l[T2] with 
effective anisotropic properties and strong coupling of the SPP modes between neighbouring layers. Otherwise 
the loss may be compensated by the use of media with gain [T31II1] . Notably, according to the the effective 
medium model, a layered lens is approximated with a uniaxially anisotropic uniform slab [12j . Uniaxial 
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Fig. 1. Schematic of a layered metal-dielectric periodic lens with N = A periods. It is 
further assumed that A = 430nm, A = 57.5nm, N = 10, c?Ag = 0.37A, dsrTiOa = 0.63A, 
UAg = 0.04 + 2.46i, nsrTiOa = 2.674 + 0.027i. 



crystals arc useful in the engineering of spatially varying polarisation states and optical vortices |151I16| and 
are a convenient medium for coupling cylindrical beams jl7j . 

Until the introduction of the concept of a perfect flat lens with either a single layer [T1[3J|31[TH] or with 
multiple layers jl3j . it was rather unusual to regard multilayers as spatial imaging systems or linear spatial 
filters. Instead, the conventional characteristics of a multilayer include the dependence of transmission and 
reflection coefficients on the angle of incidence and polarisation, as well as the photonic band structure in case 
of periodic stacks. However, in order to describe the resolution of an imaging system consisting of a multilayer 
in a systematic way, it is convenient to refer to the theory of linear shift invariant systems (LSI, also termed 
as linear isoplanatic systems [T9l[20]) and such a description already gained considerable interest [2TM25j . 
In this paper, metal-dielectric multilayers (MDM) are regarded as LSI systems and a layered superlens is a 
one-dimensional spatial filter characterised by the point spread function (PSF). This approach may facilitate 
the application of plasmonic elements to optical signal processing which is currently capturing increasing 
research interest [26| . In the present paper we introduce a vectorial description of a linear spatial filter 
consisting of a multilayer, which enables accounting for polarisation coupling that occurs during imaging. 
Therefore we examine features specific to 2D imaging when the description of the system does not trivially 
decouple into independent TE and TM polarisations. 

2. Background 

We analyse coherent imaging through a multilayer consisting of linear and isotropic materials. Further, we 
focus on metal-dielectric multilayers. A schematic of a periodic multilayered superlens is shown in Fig. [T] 
The slab boundaries arc infinite and parallel to each another, and perpendicular to the z direction. Such an 
imaging element is an example of an LSI system which realises linear spatial filtering of the TE {Ez = 0) 
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or TM {Hz = 0) -polarised incident wavefunction. Notably, diffraction in air and propagation through a 
metallic or dielectric slab are simple examples of linear spatial filters of the incident field, while a layered 
superlens is a more sophisticated one and may be seen as a superposition of such simpler systems, however 
it is still a linear spatial filter. From an engineering perspective, multiple degrees of freedom of a multilayer 
provide some leeway to perform point spread function optimisation of the response of the spatial filter. The 
distribution of scalar electric or magnetic field components in the plane of incidence 4'(a;, j/)^^^;,,^ and an 
output plane '^{x^y)z=zauti with a layered structure in between, are related through a convolution relation 
with the two-dimensional point spread function 'H{x,y), 

"^{x, y)z=z^^, = {2n)-^n{x, y) * (1) 

Here * denotes a two-dimensional convolution and we assume that propagation is in the z direction. The 
same relation transformed to the domain of spatial harmonics provides the definition of the transfer function 
(TF) which we denote as H, 

^{k^,ky)z=z„^, = 'H{k^,ky) ■ i!{K,ky)z=z,„,, (2) 
where the hat represents the two-dimensional unitary Fourier transform, defined as 

^{x,y)z ^(271)-' I J i'{k,,ky)z e'^'^-^+^-yUk^dky. (3) 

The PSF, related to the TF through the spatial Fourier transform, provides precise information on reshaping 
of the optical signal due to diffraction combined with multiple reflections and refractions, and is direct and 
convenient to interpret. Resolution of the system depends on the bandwidth of the transmitted spatial 
spectrum and super-resolution in the near field occurs when the transfer function of the system does not 
vanish for the evanescent part of the spatial spectrum (fc^ > 1). 

Until now, we assumed that the system is scalar. This is only true when the field is purely TE or TM- 
polarised or when the transfer function is the same for the TE or TM polarisations which basically occurs 
only for simple propagation in air. There exist at least two important situations when the field symmetry 
in our system decouples Maxwell's equations into these two polarisation states. One takes place for in-plane 
imaging (e.g. in the x — z plane) of ID field distributions, when imaging of scalar fields {Hy{x), Ex{x), Ez{x)) 
and {Ey{x), Hx{x), Hz{x)) are independent and stand for the TM and TE polarisations, respectively. The 
other important situation takes place for a cylindrical beam with the angular wavenumber equal to m = 0, 
when the radial {Er{r), H^{r), Ez{r)) and angular {Hr{r), Eif,{r), Hz{r)) polarisations are also examples of 
pure TM and TE polarisations. In both situations, the symmetry of the system requires that the incident 
field is only one-dimensional with the dependence on either x or r. More complex incident field distributions 
are neither TE or TM polarised and due to the different transfer function for the TE and TM polarisations, 
are subject to a change in the polarisation state during imaging through a multilayer. 

3. 2D imaging properties of layered systems in 3D 

At a plane z = zq, a monochromatic field is characterised by two polarisation components. For in- 
stance, for the linear or circular polarisation these are Ej.(x,y)z=zo ^^'^ Ey{x,y)z=zo^ o'" Ecr^{x,y)z=zo and 
Ea. {x, y)z=zo , respectively. 

Initially, let us assume knowledge of the one-dimensional transfer function (TF) for the TE and TM 
polarisations, denoting them as i-LTsikx) and HxMikx), respectively. These TFs may be easily determined 
using the transfer matrix method from the complex transmission coefficients of propagating and evanescent 
planewaves with the wavevector components in the x — z plane. We introduce denotations for the mean and 
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difference of these two TFs, Tim = CHtm +'Ht£;)/2, and Hs = {V-tm — 'Ht_e)/2, which are more convenient 
for the description of 2D imaging. 

More generally, the incident and output fields are not scalar but vectorial and the 3D transmission of 
2D field distributions through the multilayer transforms the spatial distribution of the polarisation state. 
Obviously in a planar geometry, the description of the system decouples into TE and TM polarisations in 
a simple way. In 3D, the polarisation-dependent 2D transfer function H and point spread function Ji, take 
the form of 2 x 2 matrices responsible for the transfer and coupling of both polarisations. We will now derive 
these functions for linear and circular polarisations. In cylindrical coordinates, the following polarisation 
dependent relations can be written between the incident and output fields for linear polarisations. 
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respectively, where "Hunikr, k^p) and T-lcircikr, kip) are TFs in the matrix form which we need to determine. 
The right sides of Eqs. (j4]) and ([5]) include matrix-vector multiplications. Notably, using the Jones calculus, 
it is straightforward to modify Eqs. Q and ^ and obtain expressions for the cross-coupling between linear 
and circular polarisations as well. Now, we will express these 2D TFs by the TFs for the TE and TM 
polarisations, by rotating the coordinate system and casting the 3D situation to a planar TE or TM case. 
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where R is the rotation matrix and Hxxikr, kip) may be understood as the transfer function linking the 
x-polarised component of the input signal to the x-polarised component of the output signal, T-Lxyikr, kip) as 
the transfer function linking the x-polariscd input signal to the y-polarised component of the output signal, 
etc. The mean value of the TE and TM TFs T-Lm{kr) stands for the angularly independent term of the 
transfer function T-Lun- On the other hand, the difference between the TM and TE TFs, T-Ls{kr) appears in 
the transfer function l-Lun with an angular dependence and is responsible for polarisation coupling. This term 
vanishes when TE and TM polarisations become degenerated, for instance for diffraction in air. Conversely, 
plasmon-assisted transmission of evanescent waves is possible for TM polarisation only, and this term then 
results in strong polarisation coupling in 2D. 

A similar expression for the transfer matrix for circular polarisations is obtained from Eq. ([8|) using the 
relations between linear and circular polarisations, Ea^ = {Ex + iEy)/y/2 and Ea_ = {Ex — iEy)/^/2, 
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In order to derive the corresponding 2D point spread functions, we will use the following decomposition for 
the 2D Fourier transform of a function separable in the angular coordinate system. For g{r, ip) = gr{r) -g^pi^), 
the corresponding Fourier transform g{kr, k^) is equal to (20| . 

+00 

g{kr,k^)^ ^ {-iY'cn ■ exp{ink^) ■ Hn{grir)} , (11) 

n— — 00 

where, the expansion coefficients c„ are equal to 

27r 

gip{(p) ■ exp{~inip)dtp, (12) 



and H„ denotes the n-th order Hankcl transform, 

/•OO 

H« {gr{r)} = / gr{r) ■ Jn{krr) ■ r ■ dr. (13) 



This decomposition leads to the following expressions for the polarisation-dependent 2D point spread 
function of the layered system. 
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where H„(r) = Hq {H™(A:,)} and Hsir) = H2 {%(fc,.)}- 

Summarising, the numerical calculation of the point spread functions in 2D requires the following steps: 

1. finding the ID TE and TM transfer functions for a planar geometry using the transfer matrix method, 

2. determination of the ID Hankel transforms of the 0-th and 2-nd orders of the mean and difference of these 
two TFs, respectively. 3. Reconstruction of the point spread functions and ((T7)) including the angular 
dependence. Therefore, the majority of the numerical calculations is done in ID. While we evaluate directly 
the Hankel transforms, it should be noted that there exist fast numerical algorithms which may be applied 
instead [57]. The 0-th order Hankcl transform links the rotationally invariant part of PSF to the respective 
rotationally invariant part of TF, which in turn is equal to the average of TE and TM one-dimensional TFs. 
On the other hand, the difference of the TM and TE one-dimensional TFs is responsible for the angularly 
dependent polarisation coupling and is related to the respective part of the point spread function through 
the second order Hankel transform. 



4. Polarisation coupling for linearly or circularly polarised incident beams 

In this section we calculate the TF-matrix and PSF-matrix for two-dimensional vectorial spatial filtering 
realized with a layered lens. We assume that the lens consists of silver and strontium titanate with a filling 
fraction of silver equal to 0.37 and operates at a wavelength of A = 430nm. The structure is presented in 
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Fig. 2. a) One-dimensional TF and PSF of the multilayer for in-plane propagation: 
'HtMj'Hte ~ TF for the TM and TE polarisations; b) Radial cross-section of two- 
dimensional TF and PSFs of the multilayer for 3D propagation 'Hm ~ {Htm + 'Hte)/2, 
Hs ~ {'Htm — 'Hte)/2. ID and 2D PSFs and TF corresponding to propagation in air for 
the same distance are also shown for comparison. Solid lines and dashed lines represent the 
amplitudes and phases of the respective complex functions. 



Fig.[TJ It is a low-loss self-guiding superlens with a sub- wavelength FWHM of PSF for the TM polarisation in 
the order of A/ 10 and it has been already thoroughly investigated in terms of transmission efficiency as well 
as of resolution for in-planc imaging |23ll25j . The elementary cell consists of an Ag layer symmetrically coated 
with SrTiOa. Strontium titanate is an isotropic material with a high refractive index n = 2.674 4- 0.027i at 
A = 430nm [28]. The refractive index of silver at the same wavelength is equal to UAg = 0.04 -|- 2.46z [29] . 
Further, we assume that the multilayer consists of iV = 10 periods with a thickness of A = 57.5nm each. 

In Fig. [5] in the top row we show the phase (dashed lines) and amplitude (solid lines) of the TFs: 
T~Ltm{x)^'Hte{x), 'Hrn{kr), 'Hs{r), and the TF for propagation in air at the same distance. In the bottom row, 
the corresponding ID PSF and the radial part of the 2D PSF are also presented. Our example illustrates some 
of the general properties of the TF - for normal incidence {k^ = 0), we have 'Htm{Q) = ^T£;(0) = ?^m(0), and 
T~Ls{r) = 0. On the other hand, for evanescent harmonics SPPs are only supported for the TM polarisation 
and therefore for a sufficiently large we have iHTMikx)] >> |'HT£;(fcx)|, and in effect 2T-Ltm ~ Hm ~ Hs- 
Therefore, the bandwidth and asymptotic behaviour of Hm and H.s arc the same as of H.tm (sec the solid line 
in Fig.[2j3, top for kr/ko > 1.8). If the bandwidth of Htm extends to fc^; >> 1 and its phase is approximately 
constant, which are the necessary conditions for super-resolution, the same properties are also satisfied by 
Hm and Hs- Therefore, we expect to see super-resolution in 2D, even though 2D-imaging depends on both 
TE and TM polarisations and the TE polarisation itself does not allow for super-resolution. This observation 
is confirmed by the shape of 2D PSF presented in the same figure. We further see (the violet line in Fig. [21d, 
bottom) that the polarisation coupling determined by Hs is zero at r = and smoothly increases with kr 
reaching a maximum at r/A ~ 0.18 i.e. at (approximately) the first minimum of Hm- 

In Figs. [3] and H] we show the same 2D TF and PSF as in Fig.[2]3, plotted however in the 2D space, including 
their angular dependence according to Eqns. (|15|) and (|17|) . The polarisation-preserving and polarisation- 
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Fig. 3. Two-dimensional TF-matrix of the multilayer, a) Matrix elements T-Lxxik^, ky) and 
T^xyikx, ky) for linearly polarised light; b) Matrix elements 'Ha+a-^ {kx, ky) and T-La+a- {kx^ky) 
for circularly polarised light. Phase isolines are drawn at the distances of 7r/2 and are changed 
between solid and dashed lines every tt. 




a) b) 

Fig. 4. Two-dimensional PSF-matrix of the multilayer, a) Matrix elements T-Lxxixjn) and 
T~Lxy{x^y) for linearly polarised light; b) Matrix elements 'Ha^„^{x,y) and T-La^a^ix^y) for 
circularly polarised light; Phase isolines are drawn at the distances of 7r/2 and are changed 
between solid and dashed lines every tt. 
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Fig. 5. a) Contributions to the axial Poynting vector from both polarisations at the output 
plane, resulting from a linearly x-polarised incident point source; b) Axial component of the 
Poynting vector in the x — z and y — z cross-sections inside the structure normalised with 
respect to Pz{x = 0,y = 0, z = L). The simulations were performed using FDTD. The (red) 
lines separate areas with positive and negative direction of the power flow. 



Intensity and polarisation inside the mask 




Intensity and polarisation at the output plane of the multilayer 
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Fig. 6. Intensity of the electric field and the state of polarisation inside the mask (a) and at 
the output plane of the multilayer (b). The mask consists of a perfect conductor with two 
circular slits and is illuminated with a linearly x-polarised Gaussian beam. FDTD simulation, 
linear intensity mapping. 



Fig. 7. Time-averaged Poynting vector component {Pz) in the x — z cross-section of the 
structure (where x = 0). The (red) hues separate areas with the positive and negative 
values of and the (white) arrows show the orientation of the energy flow in the vertical 
direction. 



coupled elements of the TF and PSF matrix are drawn for both linear and circular polarisations. Notably, 
the size of the PSF is sub-wavelength and the polarisation state is preserved in the central part of the PSF 
because the second order Hankel transform vanishes at the origin and the strongest coupling to the orthogonal 
polarisation is observed at some distance from the center. For circular polarisation, the term 'Haj^a+{x,y) 
depends on T-Lmir) and is rotationally symmetric while the term Ha^a^ix^y) depends on 'Hs{r) and has 
a simple angular phase dependence (see Eqll71 and Fig. IHd). For linear polarisation, the term Hxxix^y) 
depends on both Hm{r) and Hs{r) (see EqUSl and Fig. |4^) and is not rotationally symmetric. 

We have also calculated the same elements of the PSF-matrix using body-of-revolution FDTD and obtained 
good agreement with the previous results calculated using the transfer matrix method and Hankel transforms. 
For comparison, the FDTD simulations arc shown in Fig. [5] where we plot the axial energy flux resulting 
from a linearly x-polarised incident point source. Figure [5^ illustrates the contributions from the two linear 
polarisations to the Poynting vector P^ at the output plane. These contributions may be directly matched to 
the 2D PSF presented in Fig.Hh- Figure[5}D shows the Poynting vector P^ in the x — z and y— z cross-sections 
of the structure and shows that the direction of polarisation has a major influence on the diffraction rate in 
orthogonal directions inside the structure. The spatial distribution of Pz shown in Fig. [5K with contributions 
from both polarisations is in good agreement with Fig. |4^. 

Further, we demonstrate how the 2D PSF may be used to design a simple diffractive nanoelement. Using 
predicted properties of the 2D PSF we design a pattern that could be imprinted into a metallic mask attached 
to the layered superlens. The pattern consists of two narrow radial slits and we use the PSF to determine 
the distance between the slits so that the interference between the polarisation-coupled elements of two 
images of slits is enhanced. From the 2D PSF shown in Fig. ^jp we find that for the radial separation of 
r ~ 0.18 A the maximum of the polarisation cross-coupling term 'Hm(r) coincides with a minimum of the direct 
coupling term 'Hs{r). The FDTD simulation showing the operation of the proposed element is presented in 
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Figs, ini and [71 A linearly polarised Gaussian beam is diffracted on the two circular slits in the mask. For 
the purpose of simplicity the mask is made of a perfect conductor (PEC) and the width of slits is equal to 
4 nm which corresponds to 8 grid-points of the FDTD 2D mesh. Inside the slits the wave becomes polarised 
perpendicularly to the slits. This is shown in Fig. It should be noted that inside the slits the field is 
neither TE or TM polarised and 7^ 0, if ^ 7^ 0. The polarisation at the output plane is shown in Fig. [S]3. 
In between the images of two slits there is a local maximum of intensity resulting from the interference of the 
2D PSFs with the cross-polarised term eliminated. At the external side of the images of the slits, without 
interference, the cross-polarised terms of PSF are still present and the polarisation state varies considerably 
with respect to that within the slits. At the same time, in between the slits, the phase shift of Hmir) reaches 
approximately tt with respect to its central value, resulting in the reversal of the direction of power flow 
Pz ■ This reversal is demonstrated in Fig. [7] which shows the cross-section of the multilayer along the x — z 
axes. An inverted direction of Pz may be seen along most of the propagation distance within the multilayer, 
although within the mask itself the direction of the power is clearly defined in the forward direction only. 

5. Conclusion 

Two-dimensional beams generally do not preserve their initial spatial distribution of the state of polarisation 
after passing through a layered flat lens. This is the direct result of different transfer functions corresponding 
to the TM and TE components of the spatial spectrum and the presence of both TM and TE components in 
the spatial spectrum of a 2D beam. We use the framework of the theory of linear shift-invariant systems to 
describe the imaging system and express the PSF 2 x 2-matrix for imaging of vectorial 2D wavefronts using 
the 0-th and 2-nd order Hankel transforms of planar transfer functions for the TE and TM polarisations. 

We calculate the polarisation-preserving and cross-polarisation 2D PSF matrix elements of a layered 
superlens showing that the resolution in 2D is of the same order as the resolution for in-plane imaging with 
the TM polarisation. This comes in spite of the fact that both the TE and TM polarised spatial harmonics 
take part in the 2D imaging. Moreover, cross-polarisation coupling is observed for spatially modulated beams 
with a linear or circular incident polarisation. 

We demonstrate how the 2D PSF may be used to design a simple diffractive nanoelement. The spatial 
separation of the PSF matrix elements for polarisation-preserving and polarisation-coupled imaging is used 
to design a mask that forms a wavcfront which during propagation exhibits destructive interference in the 
cross-polarisation term, assures the separation of non-diffracting radial beams originating from two slits in 
the mask and exhibits an interesting property of a backward power flow in between the two rings. This opens 
a possibility for further PSF engineering where the polarisation effects are accounted for. 
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